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■ ABSTRACT. The X-ray afterglow that is observed following large flares on magnetars can be 

' accurately fit by simple and quantitative theoretical models: The long term afterglow, lasting 

' of order weeks, can be understood as thermal radiation of a heated neutron star crust that 

is re-scattered in the magnetosphere. Short term afterglow is well fit by the cooling of a non- 
. degenerate, pair-rich layer, which gradually shrinks and releases heat to a pair-free zone above 

it. Measurements of persistent optical and infrared emission directly probe long lived currents 
in the magnetosphere which are a likely source of collective plasma emission. The superstrong 
magnetic field plays an important role in generating these various emissions, and previous 
inference of its strength in magnetars is supported by the good fits with observation. 

, 1. Introduction 

^ : Magnetars exhibit a variety of revealing properties: They generate soft gamma ray 

' bursts (SGR events) which typically last ~ 80 — 200 ms, and in extreme cases have 

super-Eddington luminosities, up to lO^LEdd- The brightest of these events reach peak 
luminosity within several milliseconds to a few seconds, and have observable tails that 
taper off within a few seconds to hundreds of seconds. This prompt burst emission is 
\ followed by short term X-ray afterglows lasting thousands of seconds, and, in some cases, 

' longer term X-ray afterglows which persist for weeks or longer. Three AXP's have been 

' discovered to emit optical and infra-red emission (HuUeman et al. 2000, 2001; Wang & 

Chakrabarty 2002). In one case, the optical emission is pulsed at the frequency observed 
Q ■ in the X-ray bandpass (Kern & Martin 2002) and two other AXPs have shown longer 

|H \ timescale variability (of order weeks) in their IR flux coincident with episodes of burst 

C/3 . activity (Kaspi et al. 2002a; Israel et al. 2002). 

' The light curve of some of the largest SGR bursts has been explained as a magnetically 

trapped pair fireball that shrinks to vanishing size (Thompson & Duncan 1995; Feroci 
et al. 2001). Usov (2001) has proposed that it is due to the cooling of a strange quark 
star surface just after heating by a GRB. Long term, transient X-ray afterglow has now 
been seen on two occasions from 1900-1-14, following large flares on Aug. 27, 1998 and 
April 18, 2001. The longest decay so far has been observed from 1627-41, lasting about 
3 years (Kouveliotou et al., in preparation). The recent bursts from AXP 2259-1-586 also 
show prolonged X-ray afterglow (Kaspi, Gavriil & Woods 2002a, b) . 

Transient afterglow from neutron star surfaces following episodic energy releases was 
suggested by Eichler and Cheng (1989). Such emission could be useful for probing the 
crust of the neutron star as well as the depth of the energy release. There can be 
uplifting of surface material following a sufficiently powerful release of energy just below 
the surface, which is then radiated outward as short term afterglow. If the total energy 
release is sufficiently powerful and deep, then, although most of the heat is sucked into 
the star, some transient afterglow may be observed for weeks or months. The heat that 
is absorbed by the star reemerges as steady emission over human timescales. 

Afterglow radiation also contains clues about the mechanism of SGR bursts: in par- 
ticular, it is sensitive to the temperature to which the magnetosphere is heated during 
a burst (Thompson and Duncan 1995) and the manner in which the rigid crust of the 
star yields to magnetic stresses (Lyubarsky, Eichler, and Thompson 2002). The short 
term afterglow, with an observed t^^'^ power law, can be explained as the cooling of 
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a pair-supported surface layer which is heated by exposure to an external fireball, and 
uplifted immediately thereafter. (Thompson, Woods, Eichler & Lyubarsky, in prepara- 
tion). Heat deposited at depth (or conducted inward) may also contribute significantly 
to the persistent X-ray emission of magnetars: for example, the energy of a giant flare is 
roughly comparable to the persistent emission of the SGR integrated over its observed 
history. 

The nonthermal infra-red and optical emission may be an important clue about re- 
laxation of the magnetosphere. Eichler, Gedalin & Lyubarsky (2002) suggested that the 
infra-red and optical emission may be generated by coherent plasma processes in the 
magnetospheres of magnetars much as coherent radio emission is generated in pulsar 
magnetospheres. 

2. Thermal relaxation of the crust and afterglows. 

A soft gamma ray burst may involve not only a rearrangement of the magnetic field 
outside the neutron star, but also motion, deformation, and attendant heating of the 
crust itself. The timescale of any afterglow so induced depends on the depth to which the 
crust's temperature is significantly raised by the SGR event. In order to study cooling of 
the magnetar crust, let us assume a deposition of thermal energy density of ~ 1 x 10^^ 
erg cm~^. This is near the maximum for which neutrino losses can be neglected, and 
it is comparable to the ratio of the flare energy (> 1 x 10^^ ergs) to the volume of the 
neutron star. Within the crust, this energy density is less than a percent of B^/Stt, but 
greater than the pre-existing thermal energy density at depths less than Zheat ^ 300 m 
(for a likely internal temperature of ~ 5 — 7 x lO^K; Thompson and Duncan 1996). If 
deposited over the entire surface and to a depth of ^ 500 m, this energy density implies 
a total energy of a few times the measured Aug. 27 afterglow energy. A characteristic 
feature of this heating mechanism is that the post-burst temperature increases outward 
in the heated layer, due to the strong crustal density stratification and inward heat 
conduction. 

In our model, which attributes the fading of the afterglow to the cooling of the 
magnetar surface, the key issue is the heat transfer below the surface. The super-strong 
magnetic field significantly affects the structure of the upper crust. The Landau energy is 
relativistic in a - 10^^ G magnetic field. El « SSj^^ MeV when B :$> Bqed = 4.4 x 10^^ 
G, and the electron Fermi energy, Ep, becomes comparable with the Landau energy only 
at a depth of ^ 100 m. At k;sscr depths electrons are one-dimensional. 

Below a depth of a few meters, the heat is transferred by degenerate electrons. We cal- 
culated the electron thermal conductivity making use of the code developed by Potekhin 
(1999). The electron thermal conductivity, k, has a prominent peak when Ep is about 
the Landau energy. At larger density, k decreases, reaches a minimum when electrons 
become effectively 3-dimensional (at z ^ 2zi) and then grows slowly, as in the non- 
magnetized case. At small densities (at z < zi), k rapidly decreases so that close to 
the surface the heat transfer is dominated by radiation. Close to the surface, k is so 
small that the heat resistance of the crust is dominated by the upper few meters. The 
outgoing thermal flux is formed within a "sensitivity strip" where the radiation thermal 
conductivities become comparable with the electron ones (Gudmundsson et al. 1983; 
Ventura & Potekhin 2001). 

We have developed a code for simulations of time-dependent, one-dimensional heat 
transfer within the crust of the magnetar. The calculated outgoing flux is plotted, as a 
function of time, in Fig. 1 together with the data points obtained by Woods et al. (2001). 
The initial temperature distributions in curves 1 through 4 correspond to uniform heat 
density, with T decreasing inward until it matched onto the initial (internal) value Tint- 
The heat density was normalized by the temperature T^ax at the bottom boundary 
of the skin zone. The remaining two curves show that the results are rather robust to 
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varying the initial conditions. A slight "knee" occurs when the temperature maximum 
passes the minimum of the electron conductivity (at a few times 10^ s for B — 10^^ G). 
Beyond this break, the light curve has a slope which is independent of B, because the 
thermal conductivity at greater depths approaches the B = value. An "ankle" can 
occur beyond 10^ s, when the temperature maximum merges with the interior region of 
almost constant temperature. 

We find that the transient X-ray light curve of SGR 1900+14 in the 40 days following 
the Aug. 27 event is consistent with the hypothesis that the SGR is a magnetar made 
of otherwise normal material. While there may be some freedom in choosing the heat 
deposition profile, the 40 day timescale is consistent with the basic physics of an outer 
crustal layer which is supported by relativistic degenerate electrons against gravity, and 
the heat capacity and conductivity increase considerably with depth. The power law 
index of the decay, though certainly inconsistent with a constant initial temperature, 
is found to be weakly sensitive to the exact initial temperature profile: on timescales 
more than a few days, the deeper layers are in any case cooled by inward conduction. 
Qualitatively, this causes all but ~ 20 percent of the heat to be sucked into the star 
and reradiated only over much longer timescales as surface X-ray emission or neutrinos. 
The resulting transient afterglow emission is ~ 1 percent of the flare energy, as observed 
(Woods et al. 2001), if the initial thermal energy density in the crust is comparable to 
the ratio of the flare energy to the volume of the neutron star. This is also consistent 
with the observation that the time integrated luminosity of the SGR is dominated by 
steady emission rather than by the decaying post-burst flux. 

3. Short Term Afterglow: Cooling of the Pair Supported Atmosphere 

The magnetically trapped fireball heats the surface layers of the crust and the absorbed 
heat is reradiated by surface photon emission after the the fireball has dissipated. The 
temperature of the fireball reaches about 1 MeV for the strong bursts; photons diffuse 
into the crust and heat up the surface layer of the depth about 10^ — 10^° g/cm^ (Thomp- 
son & Dimcan 1995). This temperature is just enough to dissociate nuclei (including 
^He). While the temperature is kept at 1 MeV by heating from above, the pressure is 
increased significantly over the initial hydrostatic pressure. This indicates that, in this 
conductive layer, the enthalpy per nucleon decreases with depth. After the hot magne- 
tospheric plasma dissipates, the heated layer immediately expands, and the atmosphere 
is supported by the thermal pressure of the electron-positron pairs. The temperature 
is now low enough (kT ^ 0.1-0.5 rrieC^) that helium quickly recombines, releasing up 
to ~ 7 MeV per nucleon. We assume here that most of the nuclear dissociation en- 
ergy is restored before radiative cooling is complete. (See Thompson et al. 2003 for a 
more detailed dicussion.) In this situation, the total extractable energy per nucleon s 
is constrained to a dynamic range of several: it is milikely to be more than tens of 
MeV, in the absence of fine tuning, since otherwise the material would be blown off the 
star. To summarize the above, e is likely to lie in the range of a few to tens of MeV 
per nucleon, and to decrease slowly with depth. We normalize the column of material 
to an effective Thomson depth tt = YeUT^/rap = {Y(.aT/mp) p{z')dz' and choose 

s{tt) = £o{tt/tt,o)~^- 

The atmosphere cools by photon diffusion. This process may be envisioned as a 
cooling wave, i.e. an entropy discontinuity at depth z*{t), which propagates inwards. 
We assume the strong magnetization provides stability against convection. A detailed 
analysis shows that the pair-rich zone is unstable to the formation of an annihilation 
front, above which (at z < z*) the medium is essentially pair free (Thompson et al. 
2003). Its specific enthalpy is well below that of the pair supported region, so that the 
flux which enters it from below mostly leaves the surface, and can therefore be taken 
to be constant at ^ < 2;*. At the front of the wave, the annihilating pairs release the 
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Fig. 1. Flux times 47r x 10^^ cm^ as a function of time. Curve 1 corresponds to Tmax = 5 x 10^ 
K, r,„t = 7 X 10* K, B = 10^^ G; curve 2 to T^ax = 5 x 10^ K, Tm = 7 x 10* K, B = 3 x 10^* 
G; curve 3 to T^ax = 5x10" K, Ti„t = 4 x 10* K, B = 3 x lO" G; and curve 4 to T^ax = 3 x lO'' 
K, Tint = 4 X 10* K, B = 3 X lO^"* G. The dotted curve 1' is for B = 10^" G and an initial 
temperature distribution T = 5 x lO" K at z < 30 m, and 5 x 10"if(z/30m)"" *' at z > 30 
m; curve 2' is for B = 3 X 10"G and T = 5 X 10" K at 2 < 100 m, and proportional to 
at greater depths until merging with the internal temperature of 7 x 10* K. Data points are 
from Woods et al. (2001). Squares are normalized to a distance of 9 kpc for SGR 1900+14 and 
triangles to 16 kpc. 



Pulsars, AXPs and SGRs observed with BeppoSAX and other Observatories 



5 



energy s per nucleon, so that this constant flux is just 



m„ dt 



(1) 



A detailed calculation shows that the temperature above the annihilation front varies 
slowly with depth, T ~ Tq [tt /tt.o)^ , with index 7 — ^ and normalization kT^) ^ ^m-eC^ 
at tt,o = 10^- Therefore the outgoing flux varies with time mainly because the column 
of matter above the cooling wave front increases with time. 

In the super strong magnetic field, the radiation energy transfer is dominated by 
the extraordinary photons; their free path is very high and the Rosseland mean scat- 
tering cross-section is wr itten as cr(T, B) = An'^ar {kT /meC^f{BQED / Bf (Silant'ev & 
Yakovlev 1980). The outgoing fiux is 



m„c dU Fq / B T\ 



ZYeG{B,T)dT. tt\BqedTqJ 



(2) 



where U{T) = (l/2)aT^ is the energy density of the extraordinary photons, and _fo = 
ocTq (rneC^/A:B)^/127r^. The energy balance equation (1) is now a differential equation 
for TT[z*{t)]. With the above scaling between e and tt, one finds TT[z*{t)]/TT,o = 
[(I - 5) {Y,/tI^) {B/Bqed? (FoaT/eo) i] VWs-*), and 
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For 0.4 < S < 0.75, the obtained dependence of the outgoing flux on time is close 
to the overall t^'^'^ dependence (with slight fluctuations in the spectral index of order 
0.1) observed from 1900-1-14 over 10^ s following the Aug 29, 1998 burst (Ibrahim et 
al. 2001, Lenters et al. 2003). The large temperature (kTbb 4 keV) of that afterglow 
at ^ 10 s following the burst points to a small radiative area (about 1 percent of the 



surface area of a neutron star) and a radiative flux 2 x 10 °(A;T(,i,/4 keV) erg cm ^ s 
This implies 



B 



-B, 



QED 



~ 10 ( 



lOMeV 



-3/4 / kT> 



bb 



4 keV 
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assuming 5 = 0.6. 



4. Coherent Emission from Magnetars 

Magnctars, it has been proposed (Thompson & Duncan 1996; Thompson, Lyutikov & 
Kulkarni, 2002), have twisted magnetic loops in their magnetospheres. Most of the time, 
the thermal scale height of their atmospheres is too low to populate the magnetosphere 
with thermal plasma. On the other hand, magnetospheric currents can easily be drawn 
out of the surface of the star from at least one of the footpoints. A modest rate of 
magnetic field dissipation {d\n B / dt ^ 1/lO^yr) yields a sufficient potential drop across 
the length of the loop to create enough plasma to short out any larger potential drop. 
The density of plasma so estimated is many orders of magnitude larger than that in 
pulsar magnetospheres. If pulsars can radiate coherently in the radio, this frequency 
being ultimately determined by the plasma frequency in the pulsar magnetosphere, then 
similar processes could occur in magnetar magnetospheres with the plasma frequency 
scaled up appropriately. 
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The total density must clearly be at least as high as the minimum to deliver the 
required current, 

j = -^|V X B| ~ ec X [2 X lO^^^cm'^] sin^ eA(t>N-s- (5) 
477 rig 

Here 9 is the magnetic polar angle and A(j)N-s is the relative twist (in radians) between 
the north and south magnetic poles. Because the above density greatly exceeds the 
co-rotation charge density, nearly equal numbers of positive and negative charges are 
required to avoid absurdly high electric fields. Positive charges may be supplied either 
by pulling ions off the surface of the star, or in situ through pair creation. In the second 
case, the particle density can exceed the above by a multiplicity factor rj (which may be 
quite large; e.g. Hibschmann & Arons 2001). 

That a two-species plasma is needed suggests that there is counterstreaming between 
the positive and negative charges. This gives rise to a broad band two stream instability. 
The excited plasma waves may be converted into outgoing electromagnetic waves (see, 
e.g., Gedalin, Gruman & Melrose 2002; Lyubarsky 2002). Regardless of the details of 
any particular counterstreaming model for the coherent emission, escaping coherent 
radiation probably has a frequency of the order of the plasma frequency in the frame of 
the outflowing plasma, which gives a frequency in the observer frame of 2a;p7^/^, where 
7 is the Lorentz factor, tOp = (47re^n/me)^/^, n the plasma density in the laboratory 
frame, and nie the rest mass of the electron. 

The characteristic Lorentz factor of the plasma may be estimated as that which will 
give rise to charges of both signs, which is a necessary condition for shorting out the 
strong electric fields that would otherwise obtain. To create a pair plasma in this situ- 
ation, one probably requires resonant scattering of thermal X-ray photons that emerge 
from the star's surface.^ In order to be resonantly scattered by a relativistic electron 
moving in its lowest Landau level, thermal photons of energy must have frequency 
of eB/m^c in the electron rest frame just as in sub-QED magnetic fields. Therefore a 
Lorentz factor of 7 ~ {B / BQED){meC^ / e) ~ 10^(10keV/e)i?i5 is needed. 

Making use of the above estimates for the plasma density and Lorentz factor, one 
finds the frequency for coherent emission, 

V ~ -7^a;p ~ 2 X lO" f sin^ 0A0Ar_s) ' Hz. (6) 

TT V -^6 / 

Near the surface, where Bis ^ 1 , this suggests emission in the near IR, optical, or even 
UV for high enough ry (Eichler, Gedalin and Lyubarsky, 2002). A twisted magnetic arch 
that protrudes from a magnetar surface could emit over a broad band, depending on 
the exact altitude of the emission. In the case where the current is supplied by electrons 
and ions, Lyutikov (2002) has suggested that coherent radio emission may be generated. 
(We note that strong emission at sub- millimeter wavelengths is plausible in this case.) 

The coherent emission of pulsars is only a small fraction of the spin-down power, but 
it can be a much higher fraction of the power in polar currents, as the latter is itself only 
a small fraction of the total. By the same token, a considerable fraction of the long term 
magnetic energy dissipation in magnetars could end up as coherent electromagnetic 
emission; 1 to 10 percent is not unreasonable. A magnetar at a distance of up to 10 
kpc could be detectable at 2.2 microns with imminent technology at a luminosity of 
10^^ erg/s ( ~ 10~^ of its persistent, pulsed X-ray flux). Such emission would almost 
ccrtaiuh- liave the period of the magnetar, and would probably be polarized. In analogy 

unless the magnetic field is extremely strong, B > 10-'-'' G, and strongly sheared, in -which ease a 
pair corona can be maintained through multiple non-resonant Compton scattering (Thompson et al. 
2002). 
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to pulsars, where the direction of polarization can swing with pulse phase, the time- 
integrated polarization would probably be less than that at any instant, but it could 
nevertheless be non-zero. 

Is the optical emission (Hulleman, van Kerkwijk & Kulkarni 2000) from the anoma- 
lous X-ray pulsar 4U0142+61 coherent? On the one hand, this emission has been re- 
ported (Kern & Martin 2002) to have the periodicity of the AXP, which suggests that 
it arises from near the magnetar. With correction for reddening, the optical emission 
is about 10^^{D/3Kpc)^ erg/s, and, because the emitting surface is so small, it is al- 
most certainly non-thermal. On the other hand, general energetic and thermodynam- 
ical considerations still allow incoherent optical emission from neutron star magneto- 
spheres at detectable levels (Eichler & Beskin 2000): A brightness temperature of up 
to 10^^(i3/G)~^/^ » IO^TqK is allowable in the case of emission by electrons, and this 
can in principle provide detectable optical emission from very small emitting areas. (The 
optical pulsar in the Crab exemplifies this.) The lower limit on the size of the emitting 
region is likely to come, in the case of electrons, from the constraints on the field strength 
set by the emission frequency and other considerations. The frequency and luminosity 
of the AXP optical/IR emission are also consistent with incoherent cyclotron radiation 
from a corona of hot ions beyond the radius of fast cyclotron cooling (about 30Rns; 
Thompson et al. 2002). 

Overall there are strong reasons to expect a substantial power in coherent plasma 
emission from the current- carrying magnetosphere of an AXP or SGR. Future tests 
that could conceivably establish coherence include a) ultra-fast photometry and po- 
larimetry, which could reveal rapid time variability (micropulsation) ; b) pulsed infrared 
emission, which could set much higher floors for the brightness temperature; and c) 
polarization-time profiles, which could possibly distinguish between different coherent 
emission mechanisms. The peak frequency of this emission provides a strong diagnostic 
of the energy and composition of the charge carriers. 

In conclusion, soft gamma ray repeaters and AXP's display a rich variety of transient, 
non-7-ray emission that may be caiised or otherwise affected by the SGR events. They 
seem to be very nicely explained by the the very same hypotheses that explained the 
SGR events themselves: that they are neutron stars with ultra-strong magnetic fields, 
that the SGR events are powered by magnetic energy release that can extend well into 
the neutron star crust, and that the magnctospheres are "twisted" and support long- 
lived currents. That the ultra-strong magnetic field ^ 10^''' G arises from quantitative 
fits to the data lends further support to the claim that the field is indeed so strong. 

This research was supported by an Adler Fellowship from the Israel Science Founda- 
tion, the Arnow Chair of Theoretical Astrophysics, the Israel Ministry of Absorption, 
and the NSERC of Canada. We acknowledge helpful conversations with C. Kouveliotou. 
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